Asymmetry of magnetic-field profiles in superconducting strips 
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We analyze the magnetic-field profiles H z (x) at the upper (lower) surface of a superconducting 
strip in an external magnetic field, with z perpendicular to the plane of the strip and x along its 
width. The external magnetic field H a is perpendicular or inclined to the plane of the strip. We 
show that an asymmetry of the profiles H z (x) appears in an oblique magnetic field H a and also in 
the case when the angular dependence of the critical current density j c in the superconductor is not 
symmetric relative to the z axis. The asymmetry of the profiles is related to the difference AH z (x) 
of the magnetic fields at the upper and lower surfaces of the strip, which we calculate. Measurement 
of this difference or, equivalently, of the asymmetry of the profiles can be used as a new tool for 
investigation of flux-line pinning in superconductors. 

PACS numbers: 74.25.Qt, 74.25.Sv 



I. INTRODUCTION 

In a recent Lette»i magnetic-field profiles at the up- 
per surface of a thin rectangular YBa2Cu307_a platelet 
placed in a perpendicular external magnetic field were in- 
vestigated by magneto-optical imaging, and the following 
interesting observation was made: When columnar de- 
fects slightly tilted to the c-axis (normal to the platelet 
surface) were introduced into the sample by heavy-ion 
irradiation, an asymmetry of the magnetic-field profiles 
relative to the central axis of the sample appeared, and 
this asymmetry nonmonotonically depended on the mag- 
nitude of the external magnetic field H a , disappearing at 
large H a . The authors of Ref. Q] explained the asymme- 
try by in-plane magnetization originating from a zigzag 
structure of vortices, i.e., from their partial alignment 
along the columnar defects. They also implied that when 
the vortices lose their interlayer coherence, the in-plane 
magnetization and hence the asymmetry disappear. On 
this basis, it was claimed^ that the asymmetry can be 
a powerful probe for the interlayer coherence in super- 
conductors. In this paper we show that the asymmetry 
of the magnetic field profiles in thin flat superconduc- 
tors may have a more general origin. It may result from 
anisotropy of flux-line pinning and needs not be due to 
the kinked structure of vortices. The asymmetry can 
occur not only in layered superconductors but also in 
three dimensional materials, and its disappearance can 
be understood without the assumption that the inter- 
layer coherence is lost. Interestingly, without any colum- 
nar defects, such an asymmetry of the magnetic field pro- 
files was recently observed in a NbaSn slab placed in an 
oblique magnetic field. 2 

In Ref. 1 1 1! we explained how to solve the critical state 
problem for thin flat three-dimensional superconductors 
with an arbitrary anisotropy of flux-line pinning. But 
these equations yield magnetic-field profiles in the criti- 
cal state only to the leading order in the small parameter 
d/w where d is the thickness of the flat superconductor 



and w is its characteristic lateral dimension. In this ap- 
proximation the magnetic-field component perpendicular 
to the flat surfaces of the sample, H z , is independent of 
the coordinate z across the thickness of the superconduc- 
tor and coincides with the appropriate magnetic field of 
an infinitely thin superconductor of the same shape (but 
with some dependence of the critical sheet current J c on 
H z ). As it will be seen below, the asymmetry of the pro- 
files is related to the difference AH Z of the fields H z at 
the upper and lower surfaces of the sample. Thus, to de- 
scribe the asymmetry of the H z profiles, it is necessary to 
consider these profiles more precisely, to the next order 
in d/w, taking into account the dependence of H z on the 
coordinate across the thickness of the superconductor. 

In this paper we obtain formulas for AH Z and for the 
asymmetry of the magnetic-field profiles in an infinitely 
long thin strip, and discuss the conditions under which 
the asymmetry can be observed. In particular, an asym- 
metry always will appear for strips in an oblique mag- 
netic field, and experimental investigation of this asym- 
metry provides new possibilities for analyzing flux-line 
pinning in superconductors. We also demonstrate that 
the experimental data of Ref. Q] can be qualitatively un- 
derstood by assuming some anisotropy of pinning in a 
three-dimensional (not layered) superconductor. 



II. MAGNETIC-FIELD PROFILES OF STRIPS 

In this paper we consider the following situation: A 
thin three dimensional superconducting strip fills the 
space |a;| < w, \y\ < oo, \z\ < d/2 with d <C w; a constant 
and homogeneous external magnetic field H a is applied 
at an angle 9q to the z axis (H ax = H a sm6o, H ay = 0, 
H az = H a cos 6>o ). For definiteness, we shall imply below 
that H ax is switched on first and then H az is applied, i.e., 
the so-called third scenario^ of switching on H a occurs 
(for the definition of the first and second scenarios see 
Appendix B). It is also assumed that surface pinning is 
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negligible, the thickness of the strip, d, exceeds the Lon- 
don penetration depth, and the lower critical field H c \ is 
sufficiently small so that we may put B = /j,qH. 

The symmetry of the problem leads to the following 
relationships: 



jy (x, , 



-3y( 



H z (x,z) = H z (- 



-z), H x (x, z) = H x (-x, -z), 



where j y (x, z) is the current density flowing at the point 
(x,z). In other words, the field at the lower surface of 
the strip, H~(x) = H z (x, — d/2), can be expressed via 
the field at its upper surface, Hf{x) = H z (x,d/2), as 
follows: 

H-(x)=H+(-x), 

and hence for the difference AH z (x) = H+(x) — H~(x) 
of the fields at the upper and lower surfaces we obtain 
the formula 



AH z (x)=H+(x)-H+(-x), 



(1) 



which connects this AH z (x) and the asymmetry of the 
magnetic-field profile at the upper surface of the strip. 

As mentioned in the Introduction, to leading order in 
d/w the critical state problem for such a strip can be 
reduced to the critical state problem for the infinitely thin 
strip with some dependence of the critical sheet current 
J c on H z where the sheet current is the current density 
integrated over the thickness of the strip, 

J = I j y (x,z)dz, 

J -d/2 

and J c is its critical value. The dependence J C (H Z ) re- 
sults from both a dependence of the critical current den- 
sity j c on the absolute value of the local magnetic induc- 
tion /Uo|i?| and an out-of-plane anisotropy of j c , i.e., a de- 
pendence of j c on the angle 9 between the local direction 
of H and the z axis. Since both \H\ and 6 change with z 
in strips of finite thickness, this means that J c ^ j c d 7 and 
a dependence of J c on H z appears. The function J C {H Z ) 
can be found from the equation^*** 



d = 



HZ 



dh 



H~ jc{h,H z ) 



(2) 



where the critical current density j c (H x , H z ) may have 
arbitrary dependence on the local H x and H z ; H x = 
H ax - 0.5J C (H Z ), and H+ = H ax + 0.5J C (H Z ) are the 
x component of the magnetic field at the lower and the 
upper surfaces of the strip. The function J C (H Z ) found 
from Eq. @ generally depends on the parameter H ax , 
and only within the Bean model when j c is independent 
of H, equation J3J) yields J c = j c d for any H ax . 



Formula J2J) is valid to the leading order in d/w since 
we neglected the term dH z (x,z)/dx in the equation 
rotH = j and used the expression 



dH x (x,z) 
dz 



j y (x,z) 



(3) 



with j y (x, z) = j c (H x , H z ) to derive formula (0). In this 
approximation H z is independent of z inside the strip. 
Thus AH z (x) = 0, and so the H z profile is always sym- 
metric. This profile is given by the Biot-Savart law for 
the infinitely thin strip 



H z {x)=H a 



1 

2^ 



' J{t) dt 



t 



(4) 



Here the sheet current J(x) follows from the critical state 
equations for this strip: 



J( x ) = ~j-J c [H z {x)] 



if a < \x\ < w, and 



H z {x) = 



(5) 



(6) 



when |x| < a. The points x — ±a give the position of the 
flux front in the infinitely thin strip. The integral equa- 
tions (@J - I© can be solved by either a static iterative 
method, or more conveniently by a dynamic method^i 
An analysis of H z allowing for terms of the order of 
d/w is presented in Appendix A. In this case AH Z ^ 0, 
and it is given by 



r d/2 



d r 

AH z (x) = — zj y (x,z)dz. 
ax J - d /2 



(7) 



This expression can be also derived from the following 
simple considerations: Using divH = 0, we write 



AH z{ x)=r dH ^ z) dz 



d/2 



d ' 2 dH x {x,z) 



d/2 



Ox 



dz. (8) 



From Eq. it follows that 



H x (x,z) = H x (x,-d/2) + 



-d/2 



j y {x, z')dz' . 



Inserting this expression into formula (JSJ, interchang- 
ing the sequence of the integrations, and using 
H x (x 1 —d/2) = H ax - Q.bJ(x), we find formula (7J|. It 
follows from this formula that an asymmetry of the H z - 
profile can appear only if the distribution of the current 
density j y across the thickness of the strip is asymmetric 
about the middle plane of the strip, z — 0, and if this 
distribution changes with x (the latter condition was not 
obtained in Ref. 
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III. CONDITIONS FOR ASYMMETRY 

If the angular dependence of the critical current den- 
sity is symmetric relative to the z axis, j c {—H x ,H z ) = 
j c (H x , H z ), an asymmetry of the current distribution 
across the thickness of the strip can occur only in an 
oblique applied magnetic field which breaks the relation 
H x (x,~z) = —H x (x,z). Besides this, asymmetry of the 
distribution can appear for asymmetric angular depen- 
dence of j c , j c {~H x ,H z ) ^ j c (H x ,H z ), even when the 
external magnetic field is applied along the z axis. In 
this section we consider the case of an oblique magnetic 
field. 



A. Region where H z 







For a superconducting strip in an oblique magnetic 
field, it was shown recently^ that in the region of the 
strip, |a;| < a, where a flux-free core occurs, i.e., where 
H z pa 0, the distribution of the current across the thick- 
ness of the sample is highly asymmetric even for super- 
conductors with H- independent j c (the Bean model). Us- 
ing this distribution (which depends on how the mag- 
netic field is switched on) and Eq. J7|), one can calculate 
AH z (x) in this region of the strip. Note that the results 
of such calculations may be also applied to anisotropic 
strips with j c = j c (9) since in this region of the strip the 
flux lines practically lie in the x-y plane, j c is indepen- 
dent of the coordinates, j c « j c (7r/2), and the results for 
the flux-free core obtained within the Bean model remain 
applicable to this anisotropic case. In Fig. \T]we present 
this AH z (x) for the anisotropic strip in the case of the 
third scenario^ of switching on H a when the field H ax is 
switched on first and then the component H az is applied, 
see Appendix B. Thus, a nonzero AH Z in this region of 
the sample reflects the asymmetry of the flux-free core, 
and this AH Z differs from zero in an oblique magnetic 
field for any superconductor. 

The quantity AH z (x) in the region \x\ < a depends on 
how the external magnetic field is switched on. In partic- 
ular, when H ax is applied before H az (the third scenario) 
and j c d/2 < H ax , AH z (x) is described by Eq. (B9). 
On the other hand, for the same applied field but with 
H ax and H az switched on simultaneously (the so-called 
first scenario^), we obtain from formulas of Ref. @ at 
|a;| < asin#o 



AH z {x) 



J{x) d.J(x) 
2j c dx 



(9) 



Comparison of Eqs. (B9) and 10 demonstrates that mea- 
surements of the asymmetry of the H z profiles at the up- 
per surface of the strip enable one to investigate subtle 
differences between critical states generated by different 
scenarios of switching on H a (to the leading order in the 
small parameter d/w one has H z — at |x| < a for any 
scenario). Note also that formulas of type (B9) [or 




FIG. 1: Asymmetry of the magnetic field profile AH z (x) at 
the upper surface of a thin strip with d = O.liu and jc(9) given 
by Eq. JHJ with p = 1, q = 4, see Fig. 2. Here H az = 0.4, 
Hax = 0.5; the magnetic fields are measured in units of j c (0)d. 
The dash-dot line shows 3AH z (x) calculated with Eq. 114H 
and formulas of Appendix B. The discontinuities of this line 
seen at x — 0, x = ±o, and x — ±w are due to inapplicability 
of our approximation in the vicinity of these points. For com- 
parison, the solid line with dots shows 3AH z (x) calculated di- 
rectly by solving the two-dimensional critical state problem— 
(creep exponent a — 200 S> 1; London depth A = O.ld). The 
dashed line and the solid line without dots give the profiles 
H z (x,0) and H z (x,d/2), obtained directly from this solution 
of the two-dimensional problem. 



permit one to find the sheet current J(x) in the region 
\x\ < a from AH Z . 



B. Region where H z 7^ 



Consider now the region of the strip penetrated by 
H Zl a < \x\ < w, i.e., the region where the component 
H z differs from zero. In this case it is useful to represent 
Eq. J7J) in another form, using the solution H x (z, H z {x)) 
of Eq. J2J . This solution is implicitly given by 



z + (d/2) =± 



Hx 



dh 



Hf 3c{h,H z y 



(10) 



where H~ and H x are the same as in Eq. Us- 
ing Eq. (|10f) . we can transform the integration over z 
in Eq. Q into an integration over H x . After a simple 
manipulation, we obtain 



AH *^ = -Tx 




h dh 



Uh,H g ) 



(11) 



It is clear from this formula that AH Z depends on x via 
the function H z (x), i.e., AH Z has the form: AH Z = d ■ 
(dff z /dx)F(H z ) where F is some dimensionless function 
ofH z . 
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Within the Bean model when j c is independent of H 
and J c — j c d, equation l(TT|) yields AH z (x) — 0, and 
hence the asymmetry of the profiles H z (x) in the fully 
penetrated region of the strip is absent for any H ax . The 
asymmetry appears only if there is a dependence of j c on 
|H| or if there is an anisotropy of pinning (or it may result 
from both reasons). If j c {H x , H z ) = j c (—H x , H z ), i.e., if 
the angular dependence of the critical current density is 
symmetric relative to the z axis, it follows from formula 
JI1J that AH z (x) — at H ax = 0. In other words, 
for such j c (H x , H z ) the asymmetry of the magnetic-field 
profiles can appear only in an oblique magnetic field. At 
small H ax , Eq. ifTTJl yields 



AH z (x) 



dx 



J C (H Z 



j c (H x ,H z 



(12) 



where H. x — J c (H z )/2, and J C (H Z ) is the dependence 
of the critical sheet current on H z at H ax = 0. The 
integration of this formula leads to the relationship 



V H ax 



— const 



Jc(H z ) 
jc(H x ,H z ) ' 



(13) 



which, in principle, enables one to reconstruct the func- 
tion j c (H x , H z ) (up to a constant) if the function AH z {x) 
at small H ax and the functions H z (x) and J C (H Z ) at 
H ax — are known. Here H x = J c (H z )/2, and H z = 
H z {xq). The integration in Eq. (|13f) is carried out over 
the region where H z ^ 0. 

It was shown recently^ that in the case of H-dependent 
j c the magnetic field profiles at the upper surface of the 
strip generally depend on the scenario of switching on 
the oblique magnetic field. In Ref. the profiles were 
analyzed to the leading order in d/w, and so they were 
always symmetric in x, H z {—x,d/2) = H z (x, d/2). For- 
mula l|ll|) , which describes the antisymmetric part of the 
profiles (this part appears to the next order in d/w), 
has been derived under the assumption that the sign of 
j y remains unchanged across the thickness of the strip. 
This assumption is indeed valid for the third scenario 
discussed in this paper. However, there exist scenarios, 
see, e.g., Ref. 0, when j y (x,z) changes its sign at some 
boundary z — z c {x). It is this boundary that causes the 
difference between the profiles for the different scenarios. 
Of course, the existence of this boundary also implies a 
modification of formula l|ll|) . Thus, we expect that if 
for some scenarios the symmetric parts of the H z profiles 
differ, their antisymmetric parts have to differ, too. 

The magnetic-field profiles at the upper surface of the 
stri p are ob tained either by magneto-optics, see, e.g., 
Ref. Midi llll or using Hall-sensor arrays »•* These pro- 
files enable one to find the sheet-current distribution 
J(x) i 14 i 15 i 16 i 17 i 18 i 19 i 20 and then the dependence J C {H Z )M. 
In Ref. we discussed a way how to determine j c from 
Jc(H z ). In this context, measurements of AH z {x) in 
oblique magnetic fields can provide additional informa- 
tion on the flux-line pinning in superconductors when the 
critical current density j c depends on the magnitude or 
direction of the magnetic field. 



To demonstrate this, we compare AH z (x) for two types 
of pinning: For the first type the critical current density 
depends only on the combination \H\ cos 9 = H z where 9 
is the angle between the local direction of H and the z 
axis; for the second type j c is a function of 9 only. The 
first situation occurs for the case of weak collective pin- 
ning by point defects in the small bundle pinning regime 
when the scaling approach is valid^i In this case j c de- 
pends on the combination |_ff|(cos 2 9 + e 2 sin 2 0) 1 / 2 which 
practically coincides with \ H\ cos 9 — H z if the anisotropy 
parameter e is small. Then, equation (|llfl gives 



AH z (x) = -H a 



J C (H Z ) 
dx \j c {H z ) 



0, 



i.e., the asymmetry is absent in the region |x| > a for this 
type of pinning. Here we have taken into account that in 
this situation Eq. J5J reduces to J C (H Z ) = dj c (H z ). We 
are coming now to an analysis of AH z {x) for the second 
type of pinning. 



C. j c depends only on the direction of H 

Let us consider more closely the case when j c depends 
only on 9, j c = j c (ff), where 9 is the angle between the lo- 
cal direction of H and the z axis. In other words, we shall 
analyze the situation when the dependence of j c on \H\ 
is negligible. This approximation can be justified for not 
too thick samples with anisotropic pinning^ Here we also 
assume the symmetry j c {9) = j c (—9). In this case, one 
can express the dependence of the critical current density 
j c on 9 and the dependence of J c on H z at H ax ^ in 
terms of the function J C (H Z ) at H ax = 0^ see Appendix 
C. The quantity AH z (x), Eq. 111(1 . is also expressible in 
terms of J c (H z ,0), the sheet current J C {H Z ) at H ax = 0, 



AH z {x) 



.dH, 



dx dH z 



Hi 



J c (t, 0)dt 



4t 3 



(14) 



Here we have used the parametric representation 1C1L 
and formulas (|C2|) that determine the auxiliary variables 
t + and t_. As to Eq. 112(1. it takes the form: 



A£r,(o!) dH 

= —a 



H a 



dx dH, 



Jc(H z ) 



Jc{H z )~H z {dJ c /dH z ) 



,(15) 



where J C {H Z ) is the H z dependence of J c at H ax = 0. 

We now present an example of such calculations for this 
type of pinning. Let at H ax — the following dependence 
J C (H Z ) be extracted from some experimental magneto- 
optics data: 



J C (H Z ) = j e (0)d 



1 



- p exp 



qH 

Hr-.r 



(16) 



where H cr = je(0)d/2, while j c (0) and the dimcnsionless 
p and q are positive constants. Using Eqs. I(C1(I . one can 
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FIG. 2: Angular dependence of the critical current density 
j c (0), Eq. fL7) . for p — 1, q — 4 (dashed line). The solid lines 
show the corresponding dependences of the sheet current J c 
on H z at H ax = 0, Eq. and at H ax = 0.5, 1, 2, 5. The 

dependences for nonzero H ax are obtained from Eqs. IIU2L 
llbll . The current density is measured in units of j c (0), while 
J c and H z are in units of j c (0)d. 



easily verify that the corresponding angular dependence 
of the critical current density takes the form: 

Jc(0)=Jc(0) [l+p(l + qt)exp(-qt)} , 

tan 9 = r 1 [1 + p exp(-q t)} , (17) 

where t is a curve parameter with range < t < oo. 
This dependence j c (0) is presented in Fig.|2together with 
the function J{H Z ), Eq. l(TS)). Note that for p > the 
character of this dependence j c (^) is typical of layered 
high-Tc superconductors^ i.e., j c is largest for 9 = w/2. 
The profiles H z (x) that correspond to this J C (H Z ) are 
shown in Fig. [3J In Figs. 2 and 3 we also show the 
dependences of J c on H z for H ax ^ and appropriate 
magnetic-field profiles in an oblique magnetic field. With 
the use of Eqs. ilC2(l . I|14fl - (|16l) . we calculate the quantity 
(AH z /d)(dH z /dx)~ 1 as a function of H z , Fig. H Inter- 
estingly, even at not too small H ax < 0Aj c (0)d, formula 
(11511 gives reasonable results. An example of the asym- 
metry of the 7J z -profiles at the upper surface of the strip 
is presented in Fig. 1. This asymmetry is found from the 
data of Fig. 4, the known derivative dH z (x)/dx, Fig. 3, 
and formulas of Appendix B (when \x\ < a). Note that 
the steepness of H z (x) near x — a has a pronounced ef- 
fect on the form of AH z (x) (and this steepness essentially 
depends^ on the sign of p) . The discontinuities of AH Z (x) 
at x = ±a, 0, ±u> are caused by the inapplicability of 
our formulas for AH Z there. For comparison, we also 
show AH z (x) and the profiles H z (x,0) and H z (x,d/2) 
calculated directly by solving the two-dimensional criti- 
cal state problem for a strip of finite thickness^ In this 
two-dimensional calculation the discontinuities of AH Z 
are smoothed out on scales of the order of d. 




x / a 

FIG. 3: Spatial profiles of the perpendicular field component 
H z (x) (upper plot) and of the sheet current J(x) (lower plot) 
of a thin strip with anisotropic pinning described by model 
1171 with p — 1, q = 4, see Fig. 2. The various curves corre- 
spond to increasing applied field H az = 0.4 and 1 in units of 
j c (0)d. The dotted, dashed, dot-dashed, and solid curves are 
for H ax = 0, 0.5, 1, and 2, respectively. For comparison, the 
solid curves with dots (indicating the grid) show the profiles 
for isotropic pinning (p — 0) . 



IV. STRIP WITH INCLINED DEFECTS IN 
PERPENDICULAR MAGNETIC FIELD 

In a recent Letteri the position of the so-called central 
d-line (the discontinuity line) at the upper surface of a 
thin rectangular YBa2Cu307_a platelet was measured by 
magneto-optical imaging. In this d-line the sheet current 
changes its sign and the magnetic field H z reaches an 
extremum (a minimum). When columnar defects tilted 
to the c-axis were introduced into the sample, the d-line 
shifted relative to the central axis of the platelet ji and the 
value of this shift first increased and then decreased with 
increasing perpendicular magnetic field. As mentioned 
in the Introduction, the authors of Ref. [l| explained the 
shift by the in-plane magnetization originating from a 
zigzag structure of vortices, and the disappearance of the 
magnetization and of this shift by loss of their interlayer 
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FIG. 4: The function F{H Z ) = AH x (x)(dH x /dx)- 1 d~ 1 cal- 
culated with Eqs. 11411 . 1161 for p = 1, q = 4 and various 
H ax = 0.2, 0.5, 1.1, 1.4, 2. The dots show the F(H Z ) calcu- 
lated from Eq. I15|l for H a x = 0.2. H z and H ax are measured 
in units of j c (0)d. 



coherence. It is clear that the shift of the d-line reflects 
the asymmetry of the magnetic field profiles at the upper 
surface of the sample. 

The zigzag structure of vortices occurs when the tilt 
angle A9 of the local H to the direction of the columnar 
defects is less than the so-called trapping angle 9 t ~> In 
this case a misalignment of the local H and the averaged 
direction of the kinked flux lines appears, and this mis- 
alignment is of the order of e 2 (H c i/ H)(9 t — A9) where 
H c i is the lower critical field and e is the anisotropy pa- 
rameter of the superconductor. The misalignment gener- 
ates an in-plane magnetization M x ~ e 2 H c i(0 t —A9), and 

hence one may expect that AH Z w —(d/dx) f*y 2 M x dz; 
see Eq. (JJJ. Note that if e — > (the inter layer coherence 
is lost), AH Z tends to zero. This is just the mechanism 
of the asymmetry discussed in Ref. and in this con- 
sideration the asymmetry is due to the equilibrium part 
of the in-plane magnetization. However, in our approx- 
imation, when H ^> H c i and B = /ioH, we neglect the 
misalignment, the fields of the order of H c i, and the equi- 
librium part of magnetization. In our approach we take 
into account only the nonequilibrium part of magnetiza- 
tion, which also generates an asymmetry. This asymme- 
try can occur even if H cl — ► 0. 

In Sec. Ill we have considered the case when the an- 
gular dependence of j c is symmetric about the z axis 
\je(— = jc(9)], and the asymmetry of the H z profiles is 
caused by an inclined applied magnetic field. However, 
such asymmetry can also result from the "opposite" sit- 
uation when the applied field is along the z axis while 
pinning is not symmetric about this axis. It is this situa- 
tion that occurs when columnar defects are introduced at 
an angle 9\ to the z axis. In this case j c {@) in the interval 




FIG. 5: Top: The anisotropic critical current density jc(9), 
Eq. 1181 with pi = 2, qi = 4, and tan#i = 0.5, that 
models pinning by tilted columnar defects. Bottom: Two- 
dimensional computation 6 (creep exponent a = 50 3> 1; Lon- 
don depth A = ui/30 = d/12; 120 x 28 grid points) of the criti- 
cal state in a long strip with rectangular cross section of aspect 
ratio d/2w = 0.2 and with anisotropic pinning described by 
the j c (Q) of the upper plot, exposed to a perpendicular field 
Haz = 0.2j c od. Shown are the magnetic field lines (thin lines) 
and some contour lines of the current density j(x,z)/j c o — 
-2.75, -2.25, -1.75, . . ., 2.75 (thick lines). The three inclined 
dashed lines indicate the direction 9i of the columnar pins. 
Note that the distribution of j c across the thickness of the 
strip is asymmetric; \j c (x, z) \ is maximum where the field lines 
are along the defects, and j c w j c o is nearly constant where 
the deviation from this direction is large. The line j c (x,z) = 
coincides with the central magnetic field line. 



|A0| = \6 — #i| < 0t is larger than outside this interval. 
The enhancement of j c in this interval is caused by the 
zigzag structure of vortices when a part of their length 
is trapped by strong columnar defects, while outside the 
interval pinning by the columnar defects is ineffective. 
Since flux lines are curved in the critical state, this en- 
hancement of j c leads to an asymmetry of the current- 
density distribution across the thickness of the strip, and 
thus to a nonzero AH Z . We carry out the calculation of 
the asymmetry of the H z profiles at the upper surface of 
the strip and of the shift of the d-line for the following 
dependence j c (0)'- 



jc{0) = jco [l + Pi exp[-qi(tan( 



tan< 



(18) 



which models an increased flux-line pinning by colum- 
nar defects at angles 9 near 9\ , Fig. [SJ Here p% and q\ 
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of anisotropic pinning, but there is one more reason for 
the disappearance of the shift in samples with columnar 
defects: When 9\ > 9 t and j c od/H z is less than 6% — 9t, 
the zigzag structure of vortices disappears in the sample, 
and the flux-line pinning by these defects becomes inef- 
fective. The trapping angle determining the width of the 
peak in j c decreases with increasing H as some power of 
/H if the applied field is of the order of the matching 
field (H§ is a measure of the density of the columnar 
defects) M Thus, if 9 t > 9\ at H — 0, the disappearance 
of the shift can occur at some magnetic field associated 
with This is just observed in the experiment^ Note 
that in this case the disappearance of the shift is due to 
the disappearance of the zigzag structure of vortices in 
the superconductor rather than to the loss of the inter- 
layer coherence. 

Finally, we emphasize the unresolved problem of the 
analysis presented in this section: The shift calculated for 
the experimental ratioi (d/w) ~ 0.06 is noticeably less 
than the experimentally observed shift (AH Z decreases 
with decreasing d). A variation of the parameters pi 
and q\ in Eq. I|18fl cannot change this conclusion. For 
example, although the maximum value of the function 

l^min 

(If) | increases with p 1; it tends to a limit of the 
order of (d/2) tan#i at p\ ^> 1; see Fig. 6. 



FIG. 6: Top: Asymmetric profiles of the magnetic field com- 
ponent H z (x, z) in the plane z = 0.52d slightly above the sur- 
face of the strip of Fig. 5 with anisotropic pinning by tilted 
columnar defects. Shown are the cases pi = 2 (solid lines) and 
pi = 32 (dashed lines), with qi = 4, tan# = 0.5, at applied 
field values H az = 0.02, 0.1, 0.2, 0.3, . . .. Field unit is j c0 d. 
The slightly curved, nearly vertical lines with dots are the 
locus of the minima of these H z profiles. Bottom: The posi- 
tion imin of the minimum of the profiles H z (x,0.52d) plotted 
versus the applied field H az for pi = 2 . . . 256. 



are some positive dimensionless parameters (qi ~ 9^ 2 ), 
and j c o is the current density in a sample without colum- 
nar defects (j C Q describes, e.g., pinning by point defects). 
Although formulas (JJJ and are still valid for thin 
strips with such j c {9), these formulas fail near the d-line, 
and so we carry out calculations of AH Z here, using the 
numerical solution* of the two-dimensional critical state 
problem for a strip of finite thickness. In Fig. 5 we show 
the current and magnetic-field distributions in the strip 
with pinning described by Eq. I|18|l , while in Fig. the 
H z profiles and the shift of the d-line are presented. 

In Fig. 6 the decrease of the shift with increasing H a 
can be qualitatively explained as follows: The character- 
istic angle 9 of the curved flux lines in the strip is of the 
order of j c od/H z [i.e., for most of the flux- line elements 9 
lies in the interval —(j c od/H z ) < 9 < (j c0 d/H z )]. When 
the external magnetic field increases, this angle decreases 
and tends to zero. Thus, the flux-line pinning (the cur- 
rent distribution) becomes practically uniform across the 
thickness of the sample at sufficiently large H az , and the 
shift vanishes. These considerations are valid for any type 



V. CONCLUSIONS 

We have considered H z profiles at the upper sur- 
face of a thin strip whose thickness d is much less 
than its width 2w. To the leading order in d/w, 
these profiles are symmetric in x, H z (—x,d/2) = 
H z (x,d/2). However, the analysis to the next order 
in this small parameter reveals the asymmetry of the 
profiles, H z (x,d/2) — H z (—x,d/2), which coincides with 
AH z (x) = H z (x,d/2)-H z (x,-d/2), the difference of H z 
at the upper and lower surfaces of the strip. We calculate 
AH z (x) and show that this AH Z differs from zero in an 
oblique magnetic field H a = (H ax , 0, H az ) and depends 
on the scenario of switching on this field. For definite- 
ness, we analyze in detail the so-called third scenario- 
when H ax is switched on before H az . 

In the region |x| < a where the flux- free core occurs 
[where the symmetric part of H z (x) is almost equal to 
zero], the asymmetry of the H z profiles exists even for 
the H-indcpcndent j c (the Bean model) and is due to 
the asymmetric shape of the flux-free core in the oblique 
magnetic field. Outside this region (|a;| > a) the asymme- 
try appears only if j c depends on the magnitude of the 
local magnetic induction or if there is an out-of-plane 
anisotropy of j c . The asymmetry of the magnetic field 
profiles in oblique magnetic fields was observed in Ref . 0; 
see also Fig. 3e in Ref. |23[ 

An asymmetry of the H z profiles also appears when the 
flux-line pinning is not symmetric about the normal to 
the strip plane (about the z axis). This situation occurs 
when inclined columnar defects are introduced into the 
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strip. In this context we have discussed the experimental 
data of Itaka et ali It is shown that although these data 
can be qualitatively understood from our results, there 
is a quantitative disagreement between the experimental 
and theoretical results on the d-linc shift. 
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APPENDIX A: FORMULAS FOR H z AT THE 
SURFACES OF THE STRIP 

Using the Biot-Savart law, the H z component of the 
magnetic field at the point {xq,Zq) of the strip can be 
written in the form: 



H z (xo, 20 ) 



Haz+ h 



d/2 dz r dx j»(*,*)(*-*o) . (A1) 

— d/2 J-w r 



where r 2 = (x — xq) 2 + (z 



zo 



Using the smallness 



of the ratio d/w, we now simplify this formula. Let us 
consider the integral 



Q(z,x ,z Q ) = / dxj y (x,z) 



(x ~ x ) 



1 



(x - Xq) 



(A2) 



which appears if one calculates the difference between the 
expression (|A1I) and formula (0J for the infinitely thin 
strip. The main contribution to this Q(z, x Q , z ) is deter- 
mined by the x values near xq, \x — xq\ ~ d. Since the 
current density j y (x,z) in the critical state of the strip 
changes in the x direction on a scale which considerably 
exceeds d, in the calculation of Q(z, xq, zq) we may use 
the expansion, j y (x,z) j y (xo,z) + (x — Xo)j' y (xo, z) 
where j'(xo,z) = dj y (xo, z)/8xq. Inserting this expan- 
sion into integral (|A2() . we find that 

d 2 

Q(z,x , z )K,-Tij' (x , z)\z-z \+j y (x , z)0( — ), (A3) 

and the last term in this expression may be omitted. 
Putting zq = ±d/2, we find 



±- f d/2 Q(z,x ,±^)dz = -%L± ± \&H z {x ) , (A4) 



where AH z (x) is given by formula (JJJ), and J(x) is the 
sheet current, 



J(x) = 

Eventually we arrive at 
1 



,1/2 



d/2 



j y (x,z)dz 



Hf{x) = H az 



2tt 



J(t)dt dJ(x) d ± AH z (x) 
t-x chTI 2 ' ^ ' 
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FIG. 7: The current and flux fronts in a thin strip to which 
first the in-plane magnetic field component H ax is applied 
and then the perpendicular H az (scenario 3). Shown are the 
fronts for H ax — 0, 0.15, 0.35, 0.45 and 0.5 (or larger), and 
H az = 0.25 in units of j c d. The fronts for other values of H az 



are scaled with a. The lines 
described in Appendix B. 



Zy(x), z 1 (x), Zi(x), z 2 (x) 



The first two terms in Eqs. (|A5I) yield H z (x) for the in- 
finitely thin strip, while the third and forth terms are 
corrections to this result due to the finite thickness of 
the strip. These corrections are relatively small (of the 
order of d/w). 



The above derivation of Eq. (|A4|) fails in the region of 
the strip, |x| < a, where a flux-free core occurs in the 
sample. The boundary of the core, Xf (z) [or equivalently 
Zf(x)], can be calculated from J(x) at |x| < a; see Ref.yl 
In this region, the flux lines are practically parallel to the 
surfaces of the sample, and they are sandwiched between 
the surfaces and the core. If |xo — Xf(zo)\ < d, i.e., if the 
point (xq, zq) lies near the boundary of the flux- free core, 
one cannot transform Eq. i|A2() into Eq. (|A3p . However, 
at \x\ < a we can repeat the above analysis, integrating 
over z rather than over x in expression i|A2|l . If j y (x, z) 
is almost independent of z in the region between the core 
and the surfaces of the sample, we eventually arrive at 
the same formula l|A4() for |x| < a. This formula fails 
only in the vicinity of the points x = ±a which are the 
positions of the flux front for the infinitely thin strip and 
near the points x — 0, x — ±w. 
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APPENDIX B: FLUX-FREE CORE FOR THE 
THIRD SCENARIO 

In Ref. for a thin strip in an oblique magnetic field, 
the shapes of the flux- free core and of the lines separating 
regions with opposite directions of the critical currents 
were presented only for two scenarios of switching on the 
external magnetic field (1: at constant angle 6>o, 2: first 
H az then H ax ). Below we present the corresponding for- 
mulas for the third scenario (when H ax is applied before 
H az )- As in Refill we shall use the Bean model here, i.e., 
we assume that j c is constant. However, the formulas of 
this Appendix remain true also for the case j c = j c {9) 
where 9 is the angle between the local direction of H and 
the z axis; see below. 

Let x = ±<z be the positions of the flux front in the 
infinitely thin strip. Within the Bean model, the sheet 
current J(x) at— a < x < a — w / cosh(irH az / j c d) is^S^ 



J(x) 



2 x 
— J c arctan 

7T 



V W 2 — 
y/a? — X 



(Bl) 



In Fig.[7|we show the lines z\(x) and 2:2(2;) separating the 
regions with opposite directions of the critical currents 
and the flux- free core composed of the lines z+(x) and 
z~ (x) . Using the method of Ref. H, we find the formulas 
that describe z\(x), Z2(x), z+(x) and z~(x) at —a < 
x < 0. The appropriate formulas for a > x > can be 
obtained by the substitution: 

When H ax < j c d/4 and — x\ < x < 0, one has 



_ d 
~ 2 

z~{x) 
zx(x) 



H ax J (x) 



2jc 



H„ 



Jc 

J(x) 
4jc 



d 
2' 
d 
2' 



(B2) 
(B3) 
(B4) 



while if —a < x < —x%, the z+(x) is still given by 
Eq. IE21. but 



2 7 (x) 



J(x) H a 
2j c 



Jc 



d 
2' 



(B5) 



Here J{x) is the sheet current at the point x, and the 
point x\ follows from the condition | J(x\)\ = 4:H ax . 

When jed/i < H ax < j c d/2 and — X2 < x < 0, the 
functions z+ (x), z~(x), zi(x) are described by formulas 
(IB2|) - (IB4|I where the point xi is determined by the con- 
dition 1 3(x,2)\ = 2j c d — 4H ax . At —a < x < —X2 only 
the lines z\{x), Z2{x) exist; the line z\{x) is given by 
Eq. [|B4(1 . while 



Z2(x) 



J{x) 
4j c 



(B6) 



For high values of H ax , when j c d/2 < H ax , the flux- 
free core disappears, and one has Eq. (|B4|> for zi(x) and 
Eq. for z 2 (x) at -a < x < 0. 



Using the formulas of this Appendix and Eq. (5), it is 
easy to calculate AH z (x) in the region |a;| < a. One has 



AH z (-a <x<X!) 
AH z (~xx <x<0) 
at H ax < jcd/4, 



H ax dJ 
jc dx ' 
0.25 J + H ax dJ 



2j c 



dx 



AH z (—a < x < X2 
AH z (-x 2 <x<0) = - 



_d fd£ 
~4 \~dx~. 
0.25J + H ax dJ 



dx 



(B7) 
(B8) 

(B9) 
(B10) 



at jcd/4 < H ax < jcd/2, and Eq. (B9) at j c d/2 < H ax 
in the whole interval — a < x < 0. When < x < a, one 
can use AH z (x) = —AH z (—x). 

The formulas of this Appendix are applicable to the 
case j c = j c (9) since at \x\ < a, the flux lines are prac- 
tically parallel to the strip plane and 9 ~ tt/2. Hence, 
it is sufficient to put j c = j c (^/2) in the above formulas 
and to use J(x) obtained numerically from the appropri- 
ate solution of the critical state problem for the infinitely 
thin anisotropic strip. 



APPENDIX C: j c (0) AND J c {H z ,H al 
OF J c (H z ,0) 



IN TERMS 



When j c depends only on 9, one can reconstruct this 



dependence from the J C (H Z ) obtained at H a 







3,4 



3c (9)d = J C (H Z ) - H, 



dJ c (H z 



tan# = 



dH z 
J C (H Z 



2H 7 



(CI) 



Inserting this parametric form of j c {9) into Eq. (J2J), we 
arrive at equations determining J C {H Z ) at H ax ^ Ofi 

0.5J c (H z ,H ax ) + H ax J c (t+,0) 



H z 

\0.hJ c {H z ,H ax )-H a 



2t+ '' 
J c {t-,Q) 



2t. 



1 

H 



1 



where a is the sign of [0.5Jc{H z ,H ax )-H ax ], J c (H z ,H ax ) 
denotes J C (H Z ) at a given value of H ax , and hence J c (t, 0) 
is the sheet current at H ax = 0. At H ax = one has 
t+ = i_ = H z , and Eqs. lfU2|) reduce to J c = J C (H Z ,0) 



as it should be. From the three Eqs. (|C2(1 for the three 
unknown variables J c , t + , t- {t + and are auxiliary 
variables) one finds J C (H Z ) at H ax 7^ 0. The magnetic 
field profiles in oblique applied field are then obtained by 
inserting this effective law J C {H Z ) into the critical state 
equations for the infinitely thin strip, Eqs. (QJ - 
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